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INTRODUCTION 



ABSTRACT 

The transport of photons in steady, spherical, scattering flows is investigated. The 
moment equations are solved analytically for accretion onto a Schwarzschild black 
hole, taking into full account relativistic effects. We show that the emergent radiation 
' spectrum is a power law at high frequencies with a spectral index smaller (harder 

CNl \ spectrum) than in the non-relativistic case. Radiative transfer in an expanding enve- 

lope is also analyzed. We find that adiabatic expansion produces a drift of injected 
monochromatic photons towards lower frequencies and the formation of a power-law, 
low-energy tail with spectral index —3. 

f**«» , Key words: accretion, accretion discs - radiation mechanisms - radiative transfer. 
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. It was realized long ago (Cowsik & Lee 1982 and references therein) that the divergence of the velocity field in astrophysical 
Oh' flows can provide a very efficient mechanism to transfer energy from the fluid to particles (photons, neutrinos, cosmic rays) 
q diffusing through the medium, even in the absence of shocks. In the case of photons undergoing multiple scatterings off cold 
'_h ' electrons, this effect is germane to thermal Comptonization, with the flow velocity v playing the role of the thermal velocity, 
C/3 , and is sometimes referred to as dynamical Comptonization. In a series of papers Blandford & Payne (1981a, b; Payne & 
Blandford 1981, PB in the following) were the first to emphasize the importance of repeated scatterings in a steady, spherical 
flow of depth t 1. They have shown that monochromatic photons injected in a region where ru/c ~ 1 always gain energy, 
because of adiabatic compression, and emerge with a broad distribution which exhibits a distinctive power-law, high-energy 
tail. Under the assumption that v tx r' 5 , the spectral index depends only on (3. 

Cowsik & Lee (1982), and later Schneider & Bogdan (1989), stressed that Blandford & Payne diffusion equation for the 
photon occupation number could be regarded as a particular case of the standard cosmic-ray transport equation in which 
the diffusion coefficient, jc tx r'V 7 , does not depend on the photon energy (7 = 0) and a — f3 — 2. Starting from this result, 
Schneider & Bogdan were able to generalize PB analysis to include the transition between Thomson (7 = 0) and Klein-Nishina 
(7 = 1) scattering cross-sections. The competitive role of dynamical and thermal Comptonization in accretion flows with a 
non-zero electron temperature was studied by Colpi (1988). More recently, Mastichiadis & Kylafis (1992) investigated the 
effects of dynamical Comptonization in near-critical accretion onto a neutron star. Their approach is very similar to PB, but 
the presence of a perfectly reflecting inner boundary (either the NS surface or the magnetosphere) was taken into account. 
They have shown that, in this case, the emergent spectrum is much harder than in PB and that the spectral index depends 
both on P and on the depth at the inner boundary. 

Even when dealing with black hole accretion, all previous analyses neglected both special and general relativity. The 
presence of an event horizon was not considered and only terms up to first order in v/c were retained. In this paper we extend 
PB calculations to account for relativistic effects on radiative transfer which arise when the flow velocity approaches the speed 
of light in the vicinity of the hole horizon; as in PB, we assume that scattering is elastic in the electron rest frame. The 
motivation for this work is twofold, much like in Mastichiadis & Kylafis. First: to investigate the properties of the emergent 
spectrum in a more realistic accretion scenario, in which the optical depth near the horizon is not so large to prevent photons 
from escaping. Second: to check by means of an analytical calculation the results obtained with a numerical code recently 
developed for solving the complete transfer problem in spherical flows (Zane et al. 1996). Computed spectra show, in fact, 
a power-law, high-energy tail but the spectral index depends on the optical depth at the horizon, even for fixed (3 = —1/2, 
and it is always smaller than 2 (PB result). Here we show that advection/aberration effects in the high-speed flow near the 
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horizon, due to the finiteness of the depth there, produce a power-law tail flatter with respect to PB and enable photons to 
drift also towards energies lower than the injection energy. 

In addition, we present an analysis of dynamical Comptonization in an expanding atmosphere, using the same assumptions 
of PB. We found that the solution for the emergent flux shows specular features with respect to PB. Adiabatic expansion 
produces a drift of injected monochromatic photons to lower energies and the formation of low-energy, power-law tail. In this 
case, however, the spectral index is independent on the velocity gradient and turns out to be always equal to —3. 



2 RADIATIVE TRANSFER IN A CONVERGING FLOW 

In this and in the following sections we deal with the transfer of radiation through a scattering, steady, spherical flow, 
characterized by a power-law velocity profile v oc r 13 . Under these assumptions the rest-mass conservation yields immediately 
a density profile g oc r~ 2 ~' 3 from which it follows that the electron-scattering optical depth grj (1 + (5) oc r 

The parameter (3 is positive for outflows and it has to be (3 > — 1 for the optical depth to decrease with increasing r. Units 
in which c = G = 1 are used unless explicitly stated. 

Blandford & Payne (1981a, b) and PB restricted their analysis to converging flows with a non-relativistic bulk velocity 
and to conservative and isotropic scattering in the electron rest frame. Combining the first two moment equations, written 
in an inertial frame, they found that, in diffusion approximation, the (angle-averaged) photon occupation number n obeys a 
Fokker-Planck equation. Defining r* = 3tv and looking for separable solutions of the form 

n{v,T*) = f{T*)T*' A+ Pv-\ (1) 

the resulting second order ordinary differential equation for / reduces to a confluent hypergeometric equation (here the sign 
of (3 is the opposite with respect to PB, e.g. /3 = —1/2 for free-fall, according to the assumptions at the beginning of this 
section). As discussed by PB, the solution corresponding to a constant radiative flux at infinity and to adiabatic compression 
of photons for r — > oo is expressed in terms of the Laguerre polynomial L^ +i3 \t*). The above two conditions give rise to a 
discrete set of eigenvalues for the photon index A which are given by 

A "= 3 v + y (-=0-1.2,...) (2) 

The general solution is written as the superposition of different modes. Assuming that monochromatic photons with 
v = are injected at r* = To, it is 

oc . 
n = 

In this case, bulk motion comptonization tends to create a power law, high energy tail. Defining the spectral index as 

_ rilnL^O) 
d In v 

where L is the luminosity, PB found 

showing that the spectral slope at high frequencies is dominated by the fundamental mode n = 0. In particular, it is a = 2 
for a free-falling gas. 

The same results can be recovered using the PSTF moment formalism introduced by Thorne (1981), as shown by Nobili, 
Turolla & Zampieri (1993). Here we outline the general method, mainly to introduce same basic concepts that will be used in 
the following sections. In particular, we consider the first two PSTF moment equations with only Thomson scattering included 
in the source term. In the frame comoving with the fluid they read 
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(66) 



-Wl 



where a prime denotes the total derivative wrt lnr, y = ^Jl — r g /r/y/l — v 2 , r g is the gravitational radius, and v is taken 
positive for inward motion. As discussed by Turolla & Nobili (1988, see also Thorne, Flammang & Zytkow 1981), in diffusion 
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approximation the hierarchy of the frequency-integrated PSTF moments Wi is such that 
Wo 



Wi 

w 2 
w 3 



The same hierarchy can be assumed to hold also for frequency-dependent moments in a scattering medium. 

PB result can be reproduced in the limit of large r and small v, i.e. retaining in equations (6) only terms of order w , 
vwo and wo/t, and suppressing gravity, which is equivalent to set y = 1 and y' = 0. In particular, under such hypothesis, all 
terms containing both W2 and ws can be neglected and the moment equations become 

!S-«!S- 2 «+^ + «(™-y^H <*> 
•J£ - <•»-•. w 

where t = (1 + /3)t*. Equations (8a) and (8b) can be combined together to yield a second order, partial differential equation 
for the radiative flux 

Following PB, the solution of equation (9) can be found by separation of variables. Writing wi = t p hi(t)v~ a , it is easy 
to show that for p — 2 and for p = —/3, equation (9) becomes a confluent hypergeometric equation for hi. Actually the 
requirement of constant radiative flux at infinity is met only for p = 2, and in this case we get 

As previously discussed, the physical solution for wi can be obtained as a superposition of the Kummer functions 
M(—n, 3 + /?, t) oc l4 2+/3 ' (t), for n = 0, 1, . . ., with corresponding eigenvalues 

3(n + l) 
In this case 

oo 

»i = t 2 ^A n L^\t)v~ a - , (12) 

where the A n 's are constants to be fixed by the boundary conditions. At sufficiently large frequencies the spectrum is dominated 
by the fundamental mode; in particular, for /? = —1/2, it is again q = 2 . 



3 IMPORTANCE OF RELATIVISTIC EFFECTS 

Here we consider the effects of dynamical Comptonization in spherical accretion onto a non-rotating black hole, taking into 
full account both gravity and velocity terms in the moment equations. With reference to this particular problem, we can safely 
assume that matter is free-falling, v = (r/vg) 13 with (3 = —1/2. Note that in spherical accretion onto black holes the radiative 
flux is never going to influence the flow dynamics close to the horizon (see e.g. Gilden & Wheeler 1980; Nobili, Turolla & 
Zampieri 1991; Zampieri, Miller & Turolla 1996). In this case dynamics cancels, locally, gravity, so that y = 1, y' = 0. The 
moment equations look then "non-relativistic" in form, although v can be arbitrarily close to unity. Corrections due to large 
values of the flow velocity were not considered in previous works despite the fact that they are bounded to become important 
near the event horizon where v ~ 1. We note that the bulk of the emission in realistic accretion models is expected to come 
precisely from the region close to r g . As in the non-relativistic analysis presented in the last section, we consider the diffusion 
limit, truncating self-consistently both equations (6) to terms of order wq/t. The moments hierarchy, expressions (7), shows 
that all terms containing W3 can be always neglected in equation (6b), since they are of order wo/r 2 . In the present case, 
however, all other terms must be retained. In fact, vwi ~ wi ~ wo/t when v ~ 1 and W2 ~ wo/t, at least for v ~ 1> 1/r 
(see again expressions [7]). This implies that wi and W2 contribute to the same extent to the anisotropy of the radiation field. 
Note that under such conditions it is W2 = 4:ir(K — J/3) < 0, as already pointed out by Turolla & Nobili (1988), so that in 
high-speed, diffusive flows K may become less than J/3. Contrary to the case discussed in section 2 where 102 is negligible, 
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now the system of the first two moment equations is not closed. However, up to terms of order wo/r, the second moment 
equation does not contain moments of order higher than u>2 and provides then the required closure equation 

dw2 4 dwi 15 



dint 15 dint 



—VW2 

14 



4 2 3 dw 2 

15 5 14 o In v 



2 dwo 



15 "9 In f + 10 v W2 



0. 



(13) 



The complete system (6a), (6b) and (13) is awkward and a solution can be obtained only numerically. It is possible, nevertheless, 
to find an analytical solution if we consider the closure condition for w 2 which follows from equation (13) with only terms of 
order wq retained 



w 2 



4v / 
9 t { 



dwo 
din v 
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(14) 



With this closure, W2 is always negative provided that dwo/dlni' < 0. This implies that equation (14) is strictly valid only 
for tv> 1 (see expression [7]), that is to say below the trapping radius. Introducing the new dependent variables fo = vwo, 
fi = wi and h = vu>2, the moment equations become 



t dfo 1 df 
dt 2<91og^ 



of t dfl +9f 3 dh 



= 



(15a) 



tdfo 
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t_ 

th. 



,dfi 1 dh 
dt lOdlogi/ 



(I _ To) h 



+ t 



dh 
dt 



7 



h = o 



(156) 



h- 



dfo 



9t h \dlogu 



3/oJ=0, (15c) 

where t h is the value of t at the radius where v = 1, i.e. at r = r 9 in the case under examination. We note that for t h — > oo 
equations (15a, b) give exactly the low-velocity limit of PB (equations [8a, b]), irrespective of the value of v. This is because 
when th — > oo the scattering depth itself near the horizon must be very large, so the radiation field there is very nearly 
isotropic. Departures from isotropy, due both to the radiative flux wi ~rao/r and to the radiative shear W2 ~ (1/t — v)wo/t 
become vanishingly small, no matter how large velocity is. Under such conditions PB approach is still valid just because both 
wi and W2 become negligible in the moment equations, although they may be of the same order. 

The system (15) can be solved looking again for separable solutions of the type fi = gi(t)v~ a . After some manipulation, 
it can be transformed into a pair of decoupled, second order, ordinary differential equation for go(t) and gi(t), having the 
same structure. In particular for gi(t) it is 



where 



(77* + A) <7i =0 



(16) 



P = 60th , 

7 = -(20/3)t h [3t h -4 (3 + a)] , 

5 = 20th - 18 (2a + 3)t h - 40a (3 + a) , 

e = 30th [th - 4 (3 + a)] , 

rj = 10 (a + 2) t 2 h + (31a 2 /3 + 7a - 54) t h - 4a (a + 3) (a + 9) , 
A = (20/3)th [3th - 28 (3 + a)] . 

Equation (16) can be reduced to a hypergeometric equation upon the change of variables gi(t) = t p hi(z) and z = — (/3/7)t, 
where p is the solution of the quadratic equation yp 2 + (e — j)p — X — 0. A direct, but tedious, calculation shows that 



P+ = 2 
7 



9t ft 



(18) 



3th 



4(3 + a) ' 

in the limit t h — > oo, p_ =1/2 as in the case considered in the previous section, and p = p + will be used in the following to 
meet the requirement of constant radiative flux at infinity. Equation (16) can be now written in the form 



z(l-z) 



dz 2 



-+2p- 

7 



dh\ 

dz 



, t \ 5 n 
pip- 1) +p- - - 



hi=0 



(19) 



which is a hypergeometric equation. The general solution is expressed in terms of the hypergeometric function 2-Fi(a, b, c;z) 
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and the three parameters a, b, c (see Abramowitz & Stegun 1972, AS in the following, for notation) are given by the relations 



c = - + 2p 

7 

a + b+l = - + 2p 



ab = p(p - 1) +p 



S 7) 



Solving for a, b we obtain, after a considerable amount of algebra, 

_ 2 - a 2a (a + 3) 
~~ 2 W h 



(20a) 



It can be seen from equations (20) that, in the limit th — > oo, b diverges while a stays finite; in the same limit the 
hypergeometric equation reduces to the confluent hypergeometric equation (see e.g. Sneddon 1956). As discussed in section 
2, the relevant solution in the non-relativistic case is given by Laguerre polynomials and is recovered imposing a — —n, with 
n = 0, 1, . . .. The solution of equation (21) which reduces to PB for th — » oo is found imposing again that a is either zero or a 
negative integer (although other classes of solutions that do not match PB may exist). In this case hi is still polynomial and 
takes the form 



hi (z) = 2-Fi (— n, b, c; z) 



nl p(c- 
(c) " 



-n) 



(1-2*) , 



(21) 



it is c ~ 5/2, 6 ' 



where Pn P ' q (z) is the Jacobi polynomial and (c)„ = T(c + n)/r(c) is the Pochhammer's symbol (see again AS). For t h — > oo, 



,(c-l,b- 



*h/3, z 



■< 3t/t h = t/b and 



so that, as expected, the solution of section 2 is recovered. 

The discrete set of eigenvalues for the spectral index a n follows immediately from (20a) solving the quadratic equation 
a — —n. For each n both a positive, a+, and a negative, a~, mode is present 



-(12 + 3t h ) ± y/(12 + 3t h ) 2 + 96(n + l)t h 



(22) 



We checked that the eigenvalues of the equation for go are again given by equation (22), in agreement with the starting 
hypothesis that a is the same for all moments. 

The general solution for the spectral flux is obtained as a linear superposition of all modes 



Wl = t 



oo _ + oo 

V A+(-l) n( f^G n (b - n, c, z) (21) + V ^(_i)»i^ G „(6 - n, c, z) (21 

' (C)n VIA)/ ^-^ (C)n V^O 



(23) 



± 

n 5 



where we have expressed hi in terms of the shifted Jacobi polynomials G„. We remind that 6, c and « are all functions of a: 
although we dropped all indices to simplify the notation. The two sets of constants are fixed imposing a boundary condition 
at the injection frequency v — vq. The only boundary condition compatible with the assumption of a pure scattering flow for 
t < th is that all photons are created in an infinitely thin shell at t*. This is equivalent to ask that wi(t, vo) oc 8{t/t* — 1), as 
in PB. 

At variance with the results discussed in section 2, now the series in equation (23) can not be summed using the polynomial 
generating function because 6, c and z depend on n. The coefficients are solution of an upper triangular, infinite system 
of linear algebraic equations (see Appendix A). It can be easily shown that the two series in equation (23) do not converge for 
any value of v/vq. In fact, the general term of the first series, which is of the type /(n)(^o/^) Q " , can not be infinitesimal for 
arbitrarily small frequencies unless the series truncates, which is not the case if it must reproduce the 5-function at v — vq. 
On the other hand, the series is absolutely convergent for v > v , provided that \f(n)\ is bounded. For TV 3> 1, the series has 
a majorant oc y^,^L N (i'o/i>) which is convergent because (fo/f) dx is finite for v > vq. The same argument applies 
to the second series for v < vq, so that the solution satisfying our boundary condition is 



wi(t, v) 



t^A-(-ir { j^G n (b-n,c,z 



n=0 
oo 
2 \ ^ 



(c) n 
.(b)r, 



AZ(-l) n ^G n (b- 



n, c, 2 



•>(sr 
ar 1 



v < v : 



v>v {) . 



(24) 
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Figure 1. Emergent flux F v (in arbitrary units) for spherical accretion onto a Schwarzschild black hole; here = 20, t* = 0.9ij, and 
ctg = 1.54. A comparison with the analytical solution (dashed line) in the limit th — > oo is illustrated in the box for t* = 20, showing a 
good agreement with PB result. 



Equation (24) exhibits two striking features, not shared by its non-relativistic counterpart, which arise from the presence 
of advection/aberration terms in the moment equations. First of all, we note that according to equation (24) photons injected 
at v — vo can be shifted both to higher and lower energies by dynamical Comptonization. This is in apparent contrast with PB 
result that photons can only gain energy in scatterings with electrons in a converging flow (the adiabatic compression). PB 
statement is, however, correct up to 0(v) terms and their equation (8) is the low-velocity limit of the more general expression 
for the photon energy change along a geodesic (see e.g. Novikov & Thorne 1973) 

' '''' ' + \e + rjn j <j ij \ , (25) 



v 61 \ 3 

where n l is the unit vector along the photon trajectory and a 1 , 8 and Oij are the flow 4-acceleration, expansion and shear, 
respectively. In free-fall a 1 vanishes while it can be safely neglected in PB approximation being 0(v 2 ). The remaining two 
terms are both of order v 

' = 4- 

2T ' , (26) 
n n Oij — — — (3/U — 1) 

where /i is the cosine of the angle between the photon and the radial directions. The mean photon energy change can be 
obtained angle-averaging equation (25) over the specific intensity 

I„(ji) = too + Sfiwi + Y^M 2 - !)™2 + . . . . (27) 
Recalling the behaviour of the radiation moments in the diffusion limit, we get 

(28) 

The second term in square brackets arises because of shear and is negligible in PB approximation, being either 0(l/r ) or 
0(v/t). This implies that the mean photon energy change is always positive. However, when advection and aberration are 
taken into account (see equations [25], [26]) photons moving in a cone around the radial direction suffer an energy loss and 
the collective effect is stronger when the flow velocity approaches unity in regions of moderate optical depth. 

The second important feature concerns the slope of the power-law, high-energy tail of the spectrum. From equation (22) 
the fundamental mode is 

+ -(12 + 3t h ) + V(12 + 3t h f + 96th .... 

Q () = a > (29) 
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Figure 2. Emergent flux F v computed using our CRM code; the derived spectral index is 1.36. In this model th ~ 15 and the 
corresponding value of is 1.43 (dashed line). 

and, for large values of th, equation (29) gives 



At large enough frequencies the spectral index is dominated by the fundamental mode which, for t h > 1, sensibly deviates from 
the value predicted by the non-relativistic calculation. Despite the fact that this effect is present below the trapping radius, 
we stress that, contrary to a widespread belief, the trapping radius does not act as a one-way membrane. Photons produced 
near or below the surface tv = 1 can still escape to infinity even if both the large optical depth and the strong advection 
caused by the inward flow dramatically reduce the emergent radiative flux. Moreover, these photons, although comparatively 
few, are the more strongly comptonized and will anyway dominate the high-energy tail of the spectral distribution. Equation 
(30) shows that the emergent spectrum turns out to be flatter with respect to PB case. The two main features of our solution, 
harder spectrum and drift of photons below z/o, can be clearly seen in figure 1, where the emergent spectrum is shown for 



ft is interesting to compare the present, analytical solution with the numerical result obtained using the fully GR 
characteristic-ray code (CRM) described in Zane et al. (1996). In figure 2 we show the emergent spectrum relative to the 
"cold" solution for black hole accretion with g h = 1.42 x 1CP 5 gcm~ 3 . In this case both electron scattering and free-free 
emission/absorption are considered. At large enough frequencies scattering is the only source of opacity, so, in this limit, we 
expect our idealized analytical model to be representative of the realistic situation. The numerical model has th — 15 which 
corresponds to o?q = 1.43. This value is in excellent agreement with the derived spectral index a — 1.36. 

4 RADIATIVE TRANSFER IN AN EXPANDING ATMOSPHERE 

In this section we discuss the case of a pure scattering, expanding atmosphere with a power-law velocity profile 



where the subscript "»" refers to the base of the envelope, now v is taken positive outwards and PB approximation is used. 
The second order partial differential equation for the radiation flux is given by equation (9), upon the substitution of t with 
-t 



This equation could be integrated using the same technique discussed in section 2, looking for separable solutions wi = 
t 2 h\(t)v~ a . It can be easily checked that equation (32) yields again, upon factorization, a Kummer equation for hi(t), as 




(30) 



t h = 20. 




(31) 




(32) 
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in the converging flow case. A problem arises, however, as far as boundary conditions are concerned: in section 2 the only 
physically meaningful solution was selected asking that the flux become a constant for t — > and that adiabatic compression 
of photons hold for t — > oo. In that case the existence of these physical constraints was sufficient to fix univocally the 
mathematical solution. However, this particular issue turns out to be much more delicate in the wind problem. We preferred 
then to look for an alternative method of solution which allows for an easier handling of boundary conditions. Equation (32), 
describing diffusion of photons through a moving medium, is a Fokkcr-Planck equation and can be brought into the standard 
Fokker-Planck form 

<9 2 Oi) d[(-l - 13- t)m] _ dui 

dx 

j2„ 



(33) 



dt 2 dt 

where we have defined wi = t 2 ui and x = — 3/(2 + 0) Inf. The solution can be found by Fourier transforming equation (33) 
with respect to t, solving the equation for the Fourier transform u\ and then transforming back (see e.g. Risken 1989). The 
equation for the Fourier transform is obtained from equation (33) replacing d/dt by ik and t by id/dk, 



ik(l + ik) 



d In Mi d In Mi 



(l+(3)ik. 



(34) 



dik dx 

This is a first order PDE which can be solved by standard methods (see e.g. Sneddon 1957) once a boundary condition is 
specified. If we assume that monochromatic photons of frequency vo are injected at t = £*, the boundary condition for equation 
(34) is just iii = m? exp(— ikt„), here m? is the luminosity emitted at vq, and the corresponding solution is 



1+/3 



exp 



fii =u° {[l-exp(x-a- )]ifc + l} . 

j [1 — cxp(:r — xo)\ ik + 1 

The solution to equation (33) is given by the Fourier integral 



exp(a; — xo)ikt, 



(35) 



Ml 



2tt / 

J — t 



Mi exp(ikt) dk 



which can be evaluated analytically in terms of the modified Bessel function I q (see Appendix B) 



Ml =Mi 



v \ 3/2 /t,\ (2+/3)/2 



mm 



i - 



/ V \ 3 /(2+/3)' 



exp 



t + [y/v ) 



3/(2+13) 



u 



1 - (v/v ? n2+(i) 



2sJ{v/ Vo ynz+p)ut / (1 - (u/u ) 2/(3+l3) ) 



(36) 



The main advantage of solving equation (32) following the method outlined here is that the Fourier transform is automatically 
selecting the regular solution, because it can be computed only for functions that are L2 in ] — 00, oo[. In other words, it is 
the method of solution itself which is suited for finding only regular solutions and in doing so no extra constraint is required. 

The spectrum is shifted towards lower frequencies and it broadens at the same time, developing a power law, low-energy 
tail. The overall behaviour is similar to that of the converging flow but somehow reversed, since now photons can drift only to 
frequencies lower than uo. There is, however, a major difference in the power-law index a between the two cases since a does 
not depend on (3 for the wind solutions, as can be seen examining the spectral behaviour of equation (36) at low frequencies. 
Since I q (z) ~ (z/2) q /T(q + 1) when the argument is small, we have for the emergent luminosity 



L v oc 



(0 



exp 



(2+/3)/3 



-| -3-/3 



(37) 



1 - (^o) (2+/3)/S . 

if v <JC 1/0, which shows that a — —3 irrespective of the value of f3. The monochromatic flux at t — is shown in figure 3 for 
/3 = 1 and t* — 1; the power-law tail at low energies is clearly visible. 



5 DISCUSSION AND CONCLUSIONS 

In this paper we have reconsidered the transfer of radiation in a scattering, spherically-symmetric medium, extending PB 
analysis of converging flows to the relativistic case and investigating the effects of bulk motion comptonization in expanding 
atmospheres. In the low-velocity limit and assuming diffusion approximation, PB found that monochromatic photons injected 
at the base of the atmosphere always gain energy as they propagate outwards. The emergent spectrum exhibits an overall 
shift to higher frequencies and a power-law, high-energy tail with a spectral index depending on the velocity gradient. Under 
the same assumptions, the wind solution shows similar, although reversed, features. Adiabatic expansion now produces an 
overall drift toward lower energies and the formation of a power-law tail at low frequencies. In this case, however, the spectral 
index is independent on the velocity gradient and turns out to be always equal to —3. 

Both these analyses are correct to first order in v/c and can be thought to adequately describe situations in which 
bulk motion is non-relativistic in regions of moderate scattering depth. Relativistic corrections are, in fact, related to the 
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Figure 3. Emergent flux F„ (in arbitrary units) for an expanding envelope with (3=1 and t» = 1. At low energies the spectral index 
approaches —3 (dashed line). 

anisotropy of the radiation field and are washed out if r 3> 1. Obviously, if the flow is optically thin, repeated scatterings are 
ineffective no matter how large the velocity is. In outflowing atmospheres high velocities are expected at large radii, where the 
optical depth has dropped below unity, so that our assumption can be reasonable. On the other hand, in accretion flows onto 
compact objects the condition t 1 where v ~ 1 is likely to be met only when the accretion rate becomes hypercritical. This 
shows that a relativistic treatment of dynamical comptonization is indeed required in investigating the emission properties 
of accretion flows. For v ~ 1 the diffusion limit is not recovered simply asking that the radiative flux is proportional to the 
gradient of the energy density, since the radiative shear is as important as the flux. Relativistic corrections produces two 
main effects: first, photons are shifted toward both higher and lower frequencies by dynamical comptonization and, second, 
the spectrum at large frequencies is sensibly flatter than in the non-relativistic case. The spectral index now depends not 
only on the velocity gradient, but also on the value of the scattering depth at the horizon and goes to its non-relativistic 
limit when th tends to infinity. Despite the fact that relativistic effects are important only where tv > 1, that is to say below 
the trapping radius, their signature is still present in the emergent spectrum. In particular, the high energy tail is populated 
by the strongly comptonized photons coming just from this region. A similar effect was found by Mastichiadis & Kylafis 
(1992, see also Zampieri, Turolla, & Treves 1993) in an accretion flow onto a neutron star. In their case the formation of 
an essentially flat (a ~ 0) spectrum is due to the fact that photons experience a very large number of energetic scatterings 
before emerging to infinity, since no advection is present being the star surface a perfect reflector. Our spectrum is softer with 
respect to Mastichiadis & Kylafis just because a sizable fraction of the more boosted photons are dragged into the hole, but, 
at the same time, it is harder than PB since in the relativistic regime the mean energy gain per scattering is higher. From 
the mathematical point of view it is noteworthy that the assumption of a finite optical depth at the inner boundary (i.e. at 
the horizon in our model or at the reflecting surface in Mastichiadis & Kylafis) produces a fundamental mode which is flatter 
with respect to PB; in both cases PB result is recovered in the limit th — > oo. The possibility that scattering of photons in 
an accretion flow onto a black hole produces a power-law tail with spectral index flatter than 2 was also suggested in a very 
recent paper by Ebisawa, Titarchuk, & Chakrabarti (1996). Using a semi-qualitative analisys they found that the spectral 
index is close to 3/2 for large values of the optical depth at the horizon and discussed the possible relevance of this result in 
connection with the observed hard X-ray emission from black hole candidates in the high state. 

We note that for 1 < th < 32/9 the predicted spectral index (see equation [29]) is smaller than 1, implying a divergent 
frequency-integrated luminosity; this behaviour is not new and was already found by Schinder & Bogdan (1989) and Mas- 
tichiadis & Kylafis (1992). It simply reflects the fact that photons can gain an arbitrarily large amount of energy in collisions 
with the free-falling electrons. It should be taken into account, however, that, when hv w m e c 2 the electron recoil in the 
particle rest frame can not be neglected anymore, so for large enough energies our treatment is not valid, as discussed in 
more detail in Zampieri (1995). The decrease of the cross-section in the quantum limit makes the scattering process less 
efficient, producing a sharp cut-off in the spectral distribution. Finally, as already stressed by Blandford & Payne (1981a) 
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and Colpi (1988), thermal comptonization dominates over dynamical comptonization when v 2 ^S 12kT/m e . The spectral dis- 
tribution depends then on the relative strength of competitive processes such as heating/cooling by thermal comptonization 
and compressional heating and must be derived solving the radiative transfer equation in its complete form. 
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APPENDIX A 

Since the polynomials G n {b — n, c, z) appearing in equation (24) are not an orthogonal system, is not possible to derive an 
explicit expression for the coefficients A„. Here we show that these constants can be, in principle, obtained as the solution of 
an infinite system of linear algebraic equations. We note that the two sets A* are not independent, because the two expressions 
in (24) must match at v — vq, where 

wi(t,v ) = A5{t/U-l) (Al) 

(^4 is a constant related to the monochromatic flux injected at the inner boundary). Since z oc t/th and x = t/t*, the 
polynomials G n (b — n, c, z) can be expressed in terms of G„(3, 3, x), which form an orthogonal system, as 

n 

G n (&-n,c, 2 ) = ^G„ m G m (3,3,x). (A2) 

m=0 

The coefficients C„ m are solution of the upper triangular system of linear algebraic equations 
S^, . r -„M (m + 2)!(m + 2 + fc)! _ f n\ F(c + n)F(b + fc) / /3 \* 

2^ ( ~ l) [k (2m + 2)!(fc + 2)! ° nm ~ U ) T(b + n)F(c + k) I 7 *-U,...,n. (A3) 

m — k 

Recalling the standard expansion of the 5-function over an orthogonal set of eigenfunctions and using again G„(3,3,a;) as a 
basis, it is 



6(x-l)=J^ ^ + °^ G m (S,3,x). (A4) 

m=0 

Inserting (A2) and (A4) into (Al) and equating the coefficients of the polynomials of the same degree, we obtain 

V\_ir^C A ± - ( 2m + 3 ) ! m>0 fA5) 

1) {c)n G n m A n - t 2 m]{m + 2)] rn>0. (A5) 

n — m 

The numerical evaluation of A^ has been carried out truncating the series appearing in (A5) to a maximum order N ~ 60 
and solving the system by backsubstitution. 
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APPENDIX B 



In this Appendix we derive the expression for ui, equation (36), starting from the Fourier integral. By defining, for the sake 
of conciseness, a = (y /u Q )' i ^ 2+l3 \ the Fourier integral can be written as 



Ml 



/H 

= «! / 
2tt / 

J — t 



[(l-a)ifc + l] 1+/3 exp 



ikt 



aikt, 



(1 - a)ik + 1 



dk , 



(Bl) 



which can be transformed into an integral in the complex plane by introducing the new, complex, integration variable 
z = (1 — a)ik + 1: 

Hoo+l 



Ml = 



27ri 1 — a 



exp 



t + ai* 



1 - a 



/-t-ioo-f-J 
-ioo + 1 



1+/3 

z exp 



2 + Z 



1 - a 



dz . 



(B2) 



The integral appearing in (B2) defines the Bessel function of imaginary argument (see Prudnikov, Brychkov, & Marichev 
1986), 

J z 1+l3 exp \ T ^— : z + -^z- 1 } dz = 2th (^-) (2+/3V2 J 2+l3 [i2(atU) 1/2 /(l - a)] (B3) 



1-a 



-ioo+l 

so that finally we have 



1-a 



,o 



Ml = 



1-a \ t 



at s (2+/3)/2 



exp 



t + at» 



1-a 



(att ( 



\l/2 



1-a 



(B4) 



which is exactly equation (36). We note that, although the absolute convergence in the complex plane of the integral repre- 
sentation (B3) is proved only for at*/(a — 1) > 0, 2 + f3 < 1, direct substitution of (B4) into the Fokker-Planck equation (33) 
shows that (B4) is a solution with the only restriction a < 1. 
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